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Abstract 

We present an all-orders formula for the six-point amplitude of planar maximally 
supersymmetric M = 4 Yang-Mills theory in the leading-logarithmic approxima- 
tion of multi-Regge kinematics. In the MHV helicity configuration, our results 
agree with an integral formula of Lipatov and Prygarin through at least 14 loops. 
A differential equation linking the MHV and NMHV helicity configurations has a 
natural action in the space of functions relevant to this problem — the single- valued 
harmonic polylogarithms introduced by Brown. These functions depend on a single 
complex variable and its conjugate, w and w*, which are quadratically related to 
the original kinematic variables. We investigate the all-orders formula in the near- 
collinear limit, which is approached as \w\ — > 0. Up to power-suppressed terms, 
the resulting expansion may be organized by powers of log The leading term 
of this expansion agrees with the all-orders double-leading-logarithmic approxi- 
mation of Bartels, Lipatov, and Prygarin. The explicit form for the sub-leading 
powers of log \ w\ is given in terms of modified Bessel functions. 



1 Introduction 



In recent years, considerable progress has been made in the study of relativistic scattering ampli- 
tudes in gauge theory and gravity. A growing set of computational tools, including unitarity [1] , 
BCFW recursion pHH], BCJ duality [HIE], and symbology [8HT2]. has facilitated many impressive 
perturbative calculations at weak coupling. The AdS/CFT correspondence has provided access 
to the new, previously inaccessible frontier of strong coupling [T3]. The theory that has reaped 
the most benefit from these advances is, arguably, maximally supersymmetric M = 4 Yang- Mills 
theory, specifically in the planar limit of a large number of colors. Indeed, M = 4 super- Yang- 
Mills theory provides an excellent laboratory for the AdS/CFT correspondence, as well as for 
the structure of gauge theory amplitudes in general. 

One of the reasons for the relative simplicity of M = 4 super- Yang-Mills theory is its high de- 
gree of symmetry. The extended supersymmetry puts strong constraints on the form of scattering 
amplitudes, and it guarantees a conformal symmetry in position space. Recently, an additional 
conformal symmetry was found in the planar theory [T3HT9] . It acts on a set of dual variables, Xj, 
which are related to the external momenta kf by ki = Xj — Xj+i- At tree level, this dual conformal 
symmetry can be extended to a dual super-conformal symmetry [SU] and even combined with 
the original conformal symmetry into an infinite-dimensional Yangian symmetry [21]. At loop 
level, the dual conformal symmetry is broken by infrared divergences. According to the Wilson- 
loop/amplitude duality [13], [161 E]> these infrared divergences can be understood as ultraviolet 
divergences of particular polygonal Wilson loops. In this context, the breaking of dual conformal 
symmetry is governed by an anomalous Ward identity [I21I22]- For maximally-helicity violating 
(MHV) amplitudes, a solution to the Ward identity may be written as, 

^MHV = ^BDS x exp(iU (L1) 

where A® DS is an all- loop, all- multiplicity ansatz proposed by Bern, Dixon, and Smirnov [23], 
and R n is a dual-conformally invariant function referred to as the remainder function [24jrS5] . 

Dual conformal invariance provides a strong constraint on the form of R n . For example, 
it is impossible to construct a non-trivial dual-conformally invariant function with fewer than 
six external momenta. As a result, = R§ = 0, and, consequently, the four- and five-point 
scattering amplitudes are equal to the BDS ansatz. At six points, there are three independent 
invariant cross ratios built from distances xf- in the dual space, 

2 2 2 2 2 2 

_ X 13 X 46 _ S 12S45 _ ^24^15 _ s 23-556 _ ^35^26 _ S 34<S61 n „x 

U\ — — — ) u 2 — — q~ — j % — — 2 ~~ • U-^J 

X 1A X 36 S123S345 x 25 x 14 S 234S456 x 36 x 25 S 345S561 

Dual conformal invariance restricts Rq to be a function of these variables only, i.e. Rq = 
Rq(ui, 112,113). This function is not arbitrary since, among other conditions, it must be totally 
symmetric under permutations of the Ui and vanish in the collinear limit [24] . 

In the absence of an explicit computation, it remained a possibility that R$ — 0, despite 
the fact that all known symmetries allow for a non-zero function Rq(ui, u 2 , 113). However, a 
series of calculations have since been performed and they showed definitively that R$ ^ 0. The 
first evidence of a non-vanishing R$ came from an analysis of the multi-Regge limits of 2 — > 4 
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gluon scattering amplitudes at two loops [2B]. Numerical evidence was soon found at specific 
kinematic points [2HI2S], and an explicit calculation for general kinematics followed shortly 
thereafter [27|,|28]. Interestingly, the two-loop calculation for general kinematics was actually 
performed in a quasi- multi-Regge limit; the full kinematic dependence could then be inferred 
because this type of Regge limit does not modify the analytic dependence of the remainder 
function on the Wj. 

Even beyond the two-loop remainder function, the limit of multi-Regge kinematics (MRK) has 
received considerable attention in the context of M = 4 super- Yang Mills theory [25 |l25H4"Tj . One 
reason for this is that multi-leg scattering amplitudes become considerably simpler in MRK while 
still maintaining a non-trivial analytic structure. Taking the multi-Regge limit at six points, for 
example, essentially reduces the amplitude to a function of just two variables, w and w*, which 
are complex conjugates of each other. This latter point has proved particularly important in 
describing the relevant function space in this limit. In fact, it has been argued recently (40] that 
the function space is spanned by the set of single-valued harmonic polylogarithms (SVHPLs) 
introduced by Brown [42]. These functions will play a prominent role in the remainder of this 
article. 

The MRK limit of 2 — )• 4 scattering is characterized by the condition that the outgoing 
particles are widely separated in rapidity while having comparable transverse momenta. In 
terms of the cross ratios Ui, the limit is approached by sending one of the Ui, say u\, to unity, 
while letting the other two cross ratios vanish at the same rate that U\ — > 1, i.e. ui = x(l — u\) 
and m 3 = y{l — u\) for two fixed variables x and y. Actually, this prescription produces the 
Euclidean version of the MRK limit in which the six-point remainder function vanishes |43TT4o] . 
To reach the Minkowski version, which is relevant for 2 — > 4 scattering, ui must be analytically 
continued around the origin, u\ — > e~ 2m \ui\, before taking the limit. The remainder function 
may then be expanded around u\ = 1 and the coefficients of this expansion are functions of only 
two variables, x and y. The variables w and w* mentioned previously are related to x and y 

by [531133], 

1 _ ww* 

x ~ (i + w)(i + w*y y ~ (i + w)(i + w*) ' l ' J 

Neglecting terms that vanish like powers of 1 — Ux, the expansion of the remainder function 
may be written ajj]> 

< HV |mrk = 27ri £JV log"(l - Ul ) [g®(w,w*) + 2nihg>(w,w*)] , (1.4) 

1=2 n=0 

where the coupling constant for planar M = 4 super- Yang-Mills theory is a = g 2 N c / (8tt 2 ). 
This expansion is organized hierarchically into the leading-logarithmic approximation (LLA) 
with n = £ — 1, the next-to- leading- logarithmic approximation (NLLA) with n = i — 2, and in 
general the N fc LL terms with n = £ — k — 1. In this article, we study the leading-logarithmic 

1 Wc follow the conventions of ref. 1551. 
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approximation, for which we may rewrite eq. (jl.4j) as, 



C HV !lla = = ,f* ; 'Erfg&faw*) , (1.5) 
log(l - Ui) 

where we have identified rj — a log(l — U\) as the relevant expansion parameter. In LLA, the real 
part of Rq vanishes, so h\_ x {w^w*) is absent in eq. (jl.5p . Expressions for gf]_ x {w,w*) have been 
given in the literature for two, three [33], and recently up to ten [40] loops. 

An all-orders integral-sum representation for -R^ hv |lla was presented in ref. [33] and was 
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generalized to the NMHV helicity configuration in ref. [39]. (The MHV case was extended to 
NLLA in ref. [36].) The formula may be understood as an inverse Fourier- Mellin transform 
from a space of moments labeled by (u,n) to the space of kinematic variables (w,w*). In the 
moment space, ^Illa^ n ) assumes a simple factorized form and may be written succinctly to 
all loop orders in terms of polygamma functions. This structure is obscured in (w, w*) space, 
as the inverse Fourier-Mellin transform generates complicated combinations of polylogarithmic 
functions. Nevertheless, these complicated expressions should bear the mark of their simple 
ancestry. In this article, we expose this inherited structure by presenting an explicit all-orders 
formula for /^Illa directly in (w,w*) space. 

We do not present a proof of this formula, but we do test its validity using several non-trivial 
consistency checks. For example, our result agrees with the integral formula mentioned above 
through at least 14 loops. In ref. [39], Lipatov, Prygarin, and Schnitzer give a simple differential 
equation linking the MHV and NMHV helicity configurations, 

,,,* d p MHV| _ ,„ 9 R NMHV| / 1 a\ 

!h? 6 |lla_ dw 6 Illa ' ( } 

which is also obeyed our formula. In the near-collinear limit, we find agreement with the all-orders 
double-leading-logarithmic approximation of Bartels, Lipatov, and Prygarin [36] . 

This article is organized as follows. In Section [2} we review the aspects of multi-Regge kine- 
matics relevant to six-particle scattering and recall the integral formulas for Rs\lla in the MHV 
and NMHV helicity configurations. The construction and properties of single-valued harmonic 
polylogarithms are reviewed in Section |3j An all-orders expression for ^Illa is presented in 
terms of these functions in Section HJ After verifying several consistency conditions of this for- 
mula, we examine its near-collinear limit in Section [5j Section [6] offers some concluding remarks 
and prospects for future work. 



2 The six-point remainder function in multi-Regge kine- 
matics 

We consider the six-gluon scattering process g^g^ — > gig^g^g? where the momenta are taken to be 
outgoing and the gluons are labeled cyclically in the clockwise direction. The limit of multi-Regge 
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kinematics is denned by the condition that the produced gluons are strongly ordered in rapidity 
while having comparable transverse momenta, 

V\ > Vh > 2/4 > V2 , \j>l±\ ^ \P5±\ - \P4±\ - \P2±\ ■ (2-1) 

In the Euclidean region, this limit is equivalent to the hierarchy of scales, 

Sl2 > S345, S456 > S34, S45 , S 56 > S 23 , Sei, S 2 34, (2.2) 

which leads to the limiting behavior of the cross ratios (jl.2p . 

1 -Mi, tt 2 , m 3 ~ 0, (2.3) 
subject to the constraint that the following ratios are held fixed, 

x = = Oil) and y = = 0{\) . (2.4) 

L — Ui 1 — U\ 

Unitarity restricts the branch cuts of physical quantities like the remainder function Rq(ui, u 2 , u 3 ) 
to appear in physical channels. In terms of the cross ratios Ui, this requirement implies that all 
branch points occur when a cross ratio vanishes or approaches infinity. If we re-express the two 
real variables x and y by a single complex variable w, 

' and ^ n+ ^' ., , (2.5) 

(1 + W){1 + W*) [1 + W){1 + w*) 

then the equivalent statement in MRK is that any function of (w,w*) must be single-valued in 
the complex w plane. 

In the Euclidean region, the remainder function actually vanishes in the multi-Regge limit. To 
obtain a no n- vanishing result, we must consider a physical region in which one of the cross ratios 
acquires a phase ^2jo\. One such region corresponds to the 2 — > 4 scattering process described 
above. It can be reached by flipping the signs of s 12 and s 45 , or, in terms of the cross ratios, by 
rotating u\ around the origin, 

ux -»• e~ 2 ™ H . (2.6) 

In the course of this analytic continuation, we pick up the discontinuity across a Mandelstam 
cut [2"6"|I3"T]. The six-point remainder function can then be expanded in the form given in eq. (11.4)) . 

00 e-i 

< HV |MRK = 2vrz^^aMog' 1 (l-n 1 ) [g^(w, w*) + 2ttz h%\w, w*)] . (2.7) 

1=2 n=0 

The large logarithms log(l — U\) organize this expansion into the leading-logarithmic approxima- 
tion (LLA) with n = £ — 1, the next-to-leading-logarithmic approximation (NLLA) with n = £—2, 
and in general the the N fc LL terms with n = £ — k — 1. 
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In refs. [3"3"|I3"6~| an all- loop integral formula for -R^ hv |mrk was presented for LLA and NLLAl 



e R+mS \uRK = COS7TO; afe 



+ ^ E (-1)" - " / -^H^W^n) ^ 



(2.8) 

Here, u(u,n) is the BFKL eigenvalue and $R C g(i / , is the regularized impact factor. They may 
be expanded perturbatively, 

u(u, n) = -a (E Vtn + a + a 2 + (9(a 3 )) , 
$ Rcg (z/, n) = 1 + a $« g (z/, n) + a 2 n) + a 3 $g g (r/, n) + G(a 4 ) . 

The leading-order eigenvalue, E Utn , was given in ref. [2H] and may be written in terms of the 
digamma function ij)(z) = 4- logr(z), 

\Tl\ 1 \Tl\\ I n \ 

^ n = "2^V^ + V ' [ 1+W+ 2 J +^( v 1 -^+ VJ -2^(1). (2.10) 

In this article, we will only need the leading-order terms, but, remarkably, the higher-order 
corrections listed in (12. 9p may also be expressed in terms of the ip function and its derivatives [361 

m. 

Returning to (12.81) . the remaining functions are, 



Uab = -7k(«) log— = \ lK(a) logH 2 

8 u 2 8 

s = I Ik (a) hg(xy) = ~Jk(o) log 



\w\ 



(2.1i; 



and the cusp anomalous dimension, which is known to all orders in perturbation theory |47j . 

CO 

7if (a) = = 4 a - 4 C2 a 2 + 22 ( 4 a 3 - (^ Ce + 4 Cf) a 4 + • • • . (2.12) 

i=\ 

In addition, there is an ambiguity regarding the Riemann sheet of the exponential factor on the 
right-hand side of (12.81) . We resolve this ambiguity with the identification, 

-> e -*^(".») ( -J— \1±^L ) . (2.13) 



2 There is a difference in conventions regarding the definition of the remainder function. What we call R is 
called log(i?) in refs. (331.36 . Apart from the zeroth order term, this distinction has no effect on LLA terms. The 
first place it makes a difference is at four loops in NLLA, in the real part. 
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The in factor in the right-hand side of eq. ( 12 . 1 3[) generates the real parts hn in eq. f !2.Tj) . For 
example, at LLA and NLLA, the following relations (40] are satisfied^], 



L t-\ 



w, w 



h\l 2 (w,w ) = -y- ^(tfl.ii) ) + — r K ' g\_ 2 '(w,w ) log^p— ^ 2i ^ 

fc=2 

where 7^ = 4 from eq. (I2.12p . Making use of eq. (11.51) . we present an alternate form of these 
identities which will be useful later, 



Re (*rw) = (Wp- + t<> "* W ) < HV 

v 7 log(l — ui) \2 of] fj 16 / 



1 2 (2)/ *\ 1 / D MHV| \ 2 

+ r^71 9i (w,w)--(R 6 |lla 

log (1 - ux) 2 



LLA 

(2.15) 



The term proportional to g\ (w,w*) addresses the special case of £ = 2 in eq. ( 12. 14ft . 

In what follows, we will focus on the leading-logarithmic approximation of ( 12. 8p . which takes 
the form, 



< HV |lla = if £ J 



n=—oo 



00 

■ 



Ui 



(2.16) 



The ^-integral may be evaluated by closing the contour and summing residues^. To perform 
the resulting double sums, one may apply the summation algorithms of ref. [48], although this 
approach is computationally challenging for high loop orders. Alternatively, an ansatz for the 
result may be expanded around \w\ = and matched term-by-term to the truncated double sum. 
The latter method requires knowledge of the complete set of functions that might arise in this 
context. In ref. [40], it was argued that the single- valued harmonic polylogarithms (SVHPLs) 
completely characterize this function space, and, using these functions, eq. (I2.16P was evaluated 
through ten loops. 

So far we have only discussed the MHV helicity configuration. We now turn to the only 
other independent helicity configuration at six points, the NMHV configuration. In MRK, the 
MHV and NMHV tree amplitudes are equal [33142] • It is natural, therefore, to define an NMHV 
remainder function, analogous to eq. (11.11) . 

^ MHV !mrk = ^ DS x exp^MHv) • (2.17) 



3 Note that the sum over k in the formula for h!f}_ 2 would not have been present if we had used the convention 
for R in refs. [33l[36]. 

4 For the special case of n = 0, our prescription is to take half the residue at v = 0. 
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In ref. 139] , it was argued that the effect of changing the helicity of one of the positive-helicity 
gluongj was equivalent to changing the impact factor for that gluon by means of the following 
replacement, 

11 , 

> —^7- (2-18) 



-IV + § IV + f 



Referring to eq. (I2.16p . this replacement leads to an integral formula for -Rg MHV |LLA, 

" + °° dvw iv+n / 2 w* iu ~ nl2 



■ — I 

R 6 Illa = -— 2^ J 



n=—oo 



'l-u? aE ' 



u.n 



(2.19) 



Following refs. [3H] arid [10], we can extract a simple rational prefactor and write eq. (12.1 9p in a 
manifestly inversion-symmetric form, 

log(l — Ui) 1 + w* { \w w* J J 

for some single-valued functions f^'(w,w*). It is possible to obtain expressions for f^'(w,w*) 
directly from eq. (1 2 . 1 9 [) by means of the truncated series approach outlined above, for example. 
A simpler method is to make use of the following differential equation, which may be deduced 
by comparing the two expressions (12. 16[) and (12 . 1 9j) . 

* d MHV| _ 9 R NMHV| 

6 ' LLA_ 6 ' LLA ' ^ ' 

In principle, solving this equation requires the difficult step of fixing the constants of integration 
in such a way that single- valuedness is preserved. As discussed in ref. [40J, this step becomes 
trivial when working in the space of SVHPLs, which are the subject of the next section. 



3 Review of single-valued harmonic polylogarithms 

Harmonic polylogarithms (HPLs) (50] are a class of generalized polylogarithmic functions that 
finds frequent application in multi-loop calculations. The HPLs are functions of a single complex 
variable, z, which will be related to the kinematic variable w by z = —w. We will continue to 
use z throughout this section in order to make contact with the existing mathematical literature. 
In general, the HPLs have branch cuts that originate at z — —1, z = 0, or z = 1. In the present 
application, we will consider the restricted class of HPLs^l whose branch points are either z = 
or z = 1. To construct them, consider the set X* of all words w formed from the letters xq and 

5 Up to power-suppressed terms, helicity must be conserved along high-energy lines, so the helicity flip must 
occur on one of the lower-energy legs, 4 or 5. 

6 In the mathematical literature, these functions are sometimes referred to as multiple polylogarithms in one 
variable. With a small abuse of notation, we will continue to use the term "HPL" to refer to this restricted set of 
functions. 
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Xi, together with e, the empty worcfj]. Then, for each w E X*, define a function H w (z) which 
obeys the differential equations, 

d u / \ H w {z) d H w (z) 

—H X0W (z) = — - — and —H XlW (z) = _ , (3.1) 

subject to the following conditions, 

H e (z) = l, H xd (z) = -log n z, and limH w ^(z) = 0. (3.2) 

There is a unique family of solutions to these equations, and it defines the HPLs. For w ^ x^, 
they can be written as iterated integrals, 

H X0W {z) = f dz' and H XlW = f dz 1 ^1 . (3.3) 



The structure of the iterated integrals endows the HPLs with an important property: they form 
a shuffle algebra. The shuffle relations can be written as, 

H Wl (z) H W2 (z) = H ™( z )i ( 3 - 4 ) 

u>Gu>iiHM)2 

where W\mw2 is the set of mergers of the sequences w\ and u>2 that preserve their relative ordering. 
The shuffle algebra may be used to remove all zeros from the right of an index vector in favor of 
some explicit logarithms. For example, it is easy to obtain the following formula for HPLs with 
a single xi, 

H x n XlX rn = (^_^)| { ^ J ' (3 ' 5) 

After removing all right-most zeros, the Taylor expansions around z = are particularly simple 
and involve only a special class of harmonic numbers [50] , 



z l 



H mi ,...,m k (z) = ^—Z m2 ^ mk (l-1), m;>0, (3.6) 
i=i 

where Z mi ^ mk (n) are Euler-Zagier sums [5TJ[52] , defined recursively by 

n 1 

Z(n) = 1 and Z TOli ... >mjk (n) = ^ j^ Z m 2 ,...,m k (l - 1) • (3.7) 



i=i 



Note that the indexing of the weight vectors mi, ... , in eqs. (13.6p and ( 13. 7ft is in the collapsed 
notation in which a subscript m denotes m — 1 zeros followed by a single 1. 

The HPLs are multi- valued functions; nevertheless, it is possible to build specific combinations 
such that the branch cuts cancel and the result is single-valued. An algorithm that explicitly 



7 Context should distinguish the word w from the kinematic variable with the same name. 
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constructs these combinations was presented in ref. [32] and reviewed in ref. [ID]. Here we provide 
a very brief description. 

The SVHPLs C w (z) are generated by the series, 

C{z) = L x {z)L Y {z) = C ^ z ) w > ( 3 - 8 ) 

wex* 

where, 

Lx(z) = ^2 H w (z)w, L Y {z) = ^2 H^ w ){z)w . (3.9) 

Here ~ : X* — y X* is the operation that reverses words, <fi : Y* — y X* is the map that renames y 
to x, and Y* is the set of words in {2/0,2/1}, which are defined by the relations, 



yo = xo 

^(2/0,2/1)2/1^(2/0, = Z(xq,x 1 )~ x x 1 Z{xq,x 1 ), 



(3.10) 



where Z(x ,Xi) is a generating function of multiple zeta values, 

Z(xo,X!) = ^2 C(w)w. (3.11) 



wex* 



The ((w) are regularized by the shuffle algebra and obey ((w 7^ X\) = H w (l) and C( x i) — 0. 

Alternatively, one may formally define these functions as solutions to simple differential equa- 
tions, i.e. the jC w (z) are the unique single- valued linear combinations of functions H Wl (z)H W2 (z) 
that obey the differential equations 



|-£.U*) = — and ^-C XlW (z) = ^ J (3.12) 
oz z oz 1 — z 

subject to the conditions, 

C e (z) = 1 , £ x %{z) = ^[ lo S n N 2 and lim£ w ^ x n(z) = 0. (3.13) 

The SVHPLs also obey differential equations in z. Both sets of equations are represented nicely 
in terms of the generating function (13. 5)) , 

^W=(t + T^)^W |«*> = £W (f + f^) • (3.14) 

4 Six-point remainder function in MRK and LLA 

The SVHPLs introduced in the previous section provide a convenient basis of functions to describe 
the six-point remainder function in MRK. In ref. [3D], these functions were used to express the 
result through ten loops in LLA and through nine loops in NLLA. Here we use the SVHPLs to 
present a formula in LLA to all loop orders. 
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4.1 The all-orders formula 



Recall from the previous section that we defined X* to be the set of all words w in the letters xq 
and x\ together with the empty word e. Let C(X) be the complex vector space generated by X* 
and let C(£) be the complex vector space spanned by the SVHPLs, C w with w 6 X*. Denote 
by C(X)[[?y]] and C(£) [[rj\] the rings of formal power series in the variable t] = alog(l — u±) with 
coefficients in C(X) and C(£), respectively. There is a natural map, p, which sends words to the 
corresponding SVHPLs, 



p : C(X)[[r)]] -> C(£)liv]] 
w £,„ . 



(4.1; 



Using these ingredients, we propose the following formulas for the MHV and NMHV remainder 
functions in MRK and LLA, 



nMHVi 
-n-6 I LLA 



rMHV 



NMHV I 



log(l - Ui) 
2ivi 



LLA 



2m . P (xz* 

^^—p(x XZ NUHV ) 
logfl — Mi) 1 + w* V / 



w, w ) 



1 1 



(4.2) 
(4.3) 



where the formal power series X, ^C N ) MHV g C(X) [[rj\] are, 
X = 



1 — X\ 

fc-1 



•i'O 



Z 



MHV 



)2m-fe+l 



fe=l 

oo 



n=0 
fc-2 



m=0 



NMHV 



(A; — m — 1)! 

22m— fc+1 

(A; — m — 1)! 



3(n, m) ) ?? fc 
3(n,m) 1 ?7 fc 



(4.4) 



Here, the 3(w, m) are particular combinations of £ values of uniform weight n. They are related 
to partial Bell polynomials, and are generated by the series, 



exp 



k=l 



oo oo 



^^3(n,m)x> m 



n=0 m=0 



An explicit formula is, 



3(n,m) = ]T II 



(4.5) 



(4.6) 



fSeP(n,m) i 

where P(n, m) is the set of n-tuples of non- negative integers that sum to m, such that the product 
of ( values has weight n, 



P(n, m) = l{p u ■ ■ ■ , P n } Pi G N , ^ A = m, ^(2* + 



(4.7) 



i=l 



i=l 
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Similarly, an expression for the kth term of X can be given as, 

fc=0 \n=0 v 7 aeQ(n) J J 

where Q(n) is the set of integer compositions of n, 

Q(n) = |{ai,a 2 , • • • ,a m } a t E Z + , ^ a { , = . (4.9) 

Excluding the one-loop term in eq. (14. 2p . the arguments of the p functions factorize into the 
product of a £-free function, X, and a (^-containing function, ^( N ) MHV . The £-free function is 
simpler and its first few terms read, 

, f l \ f l 2 1 1 2 \ 2 

X = 1 + \-x + x 1 J T] + \-x + - X X! + - XxX + x 1 jr] 

zjj + ~ XqXi + i X XiX + i £ ^l + i XiXq + ~ XiX Xi + ^ £?£ + xfj if H . 

(4.10) 

The (^-containing functions are slightly more complicated. Their first few terms are, 

z mhv = 1 ir]+ 1 i;ro77 2 + 1 Xl3 * 1 f + ( A. x x xl- JCa^i Wh , 

1 1 1 /I 1 \ (4 - 11} 

Using eqs. (14. lOf) and (14. lip , one may easily extract for £ = 1, 2, 3, 4 (cf. eqs. (11.51) and (I4.2D ) . 
The one loop term vanishes, g$ = 0, and the other functions read, 

(2) l/i l/i l/i 

9i — 7 £0,1 + 7^ £1,0 + 2 £1,1 j 

(3) I/, l/i l/i l/i l/i l/i l/i 

#2 _ 77: £o,0,l + o ^0,1,0 + T M),l,l + 77: £l,0,0 + T A.,0,1 + 7 £l,l,0 + 7 £l.l,l , 

lo 8 4 lo 4 4 2 

#3^ = 7^ £0,0,0,1 + 777 £0,0,1,0 + 77; £0,0,1,1 + 777 £0,1,0,0 + 7 £0,1,0,1 + 7 £0,1,1,0 (4-12) 

SO 6Z 10 OZ O O 

1/1 1/1 1 1/1 1/1 1/1 

+ 7 £0,1,1,1 + 777 £l,0,0,0 + 77: £l,0,0,l + 7 £l,0,l,0 + 7 £l,0,l,l + 777 £l,l,0,0 
4 9o 12 8 4 lo 

I. I. I. 1 

+ 7 £l,l,0,l + ^ £l,l,l,0 + 2 £l,l,l,l — 77 C3 £1 • 
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Similarly, one may extract the first few /W (cf. eqs. ( TZ^Uj) and (Q]) ). finding /« = and, 

f(2) . 1 /• 

— 0^0,0,1 + 777^0,1,0 + 7^0,1,1) 

8 lb 4 

/ = 77^ £0,0,0,1 + g2 ^0,0,1,0 + g £0,0,1,1 + gg £0,1,0,0 + g £0,1,0,1 + £0,1,1,0 + 7/ £0,1,1,1 ? 

/jr\ 1 _ 1 1 1 1 

/ = O' ' ' ' 1 Y28 °' ' ' 1 ' ^ 32 °>°>o,i,i ^92 °' ' 1 ' ' iq °> >i>°>i 

+ ^ ^0,0,1,1,0 + g £0,0,1,1,1 + ^gg £o,l,0,0,0 + ^ £0,1,0,0,1 + 77^ £0,1,0,1,0 

+ g £0,1,0,1,1 + g77 £0,1,1,0,0 + g £0,1,1,0,1 + £0,1,1,1,0 + 7^ £0,1,1,1,1 _ C3 £0,1 • 

(4.13) 

We do not offer a proof that eqs. ( 14.21) and ( 14.31) are valid to all orders in perturbation 
theory. One may easily check that their expansions through low loop orders, as determined by 
eqs. (14.12p and ( 14. 13ft . match the known results [331110]. It is also straightforward to extend 
the above calculations to ten loops and confirm that the results are in agreement with those of 
ref. [ID]. Moreover, we have verified that the truncated series expansion of eq. (14.21) as \w\ — > 
agrees with that of eq. ( 12. 16ft through 14 loops. 

A comparison through such a high loop order is important in order to confirm the absence 
of multiple zeta values with depth larger than one (hereafter simply "MZVs"). To see why these 
MZVs should be absent, consider performing the sum of residues in eq. (12.1 bp . Transcendental 
constants can only arise from the evaluation the ip function and its derivatives at integer values. 
The latter are given in terms of rational numbers (Euler-Zagier sums) and ordinary C, values. 
Therefore, it is impossible for the series expansion of eq. (I2.1bp to contain MZVs. 

On the other hand, we would naively expect MZVs to appear in the series expansion of 
eq. (14. 2 p at 12 loops and beyond. This expectation is due to the fact that, for high weights, the 
y alphabet of eq. (13 . 1 0[) contains MZVs, and, starting at weight 12, these MZVs begin appearing 
explicitly in the definitions of the SVHPLs. In order for eq. (14. 2p to agree with eq. (12.1 bp . all 
the MZVs must conspire to cancel in the particular linear combination of SVHPLs that appears 
in ( 14. 2p . We find that this cancellation indeed occurs, at least through 14 loops. It would be 
interesting to understand the mechanism of this cancellation, but we postpone this study to 
future work. 

4.2 Consistency of the MHV and NMHV formulas 

The MHV and NMHV remainder functions are related by the differential equation ( I2.2ip . 

8 8 



w 



pMHVi _ „„ w pNMHVi (a 1A \ 

ti% |lla — U!— K 6 |lla- l 4 - i4 J 



dw* dw 
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Recalling that (w,w*) = (—z, —z), it is straightforward to use the formulas (13.1 4p to check that 
eqs. (14. 2p and (14. 3 p obey this differential equation. To see how this works, consider eq. (14. 2p . 
which we write as, 



MHVl 



2ixi 



I LLA 



log(l - Ui) 



P 



gQ{x ,xi)xQ + g l {x Q ,x 1 )x l 



(4.15) 



for some functions ^(^o^i) an d gi(x ,Xi) which can be easily read off from eq. (14. 2p . The w* 
derivative acts on SVHPLs by clipping off the last index and multiplying by 1/w* if that index 
was an Xo or by —1/(1 + w*) if it was an x\. There are also corrections due to the y alphabet at 
higher weights. Importantly, y$ = Xq, so these corrections only affect the terms with a prefactor 
1/(1 + w*). This observation allows us to write, 



w 



u pMHV 

dw* 6 



2m 



LLA 



log(l -m) 

2wi 
log(l - Ui) 



-P 



1 



w 



w* 



1 + w 

goixoiXx) -+ 



:g 1 {X ,X 1/ 

1 

1 + 1/w* 



(4.16) 



Due to the complicated expression for yi, it is difficult to obtain an explicit formula for gi(xo, x±). 
Thankfully, we may employ a symmetry argument to avoid calculating it directly. Referring 
to eq. (I2.16p . -R^ hv |lla has manifest symmetry under inversion (w,w*) O (1/w, 1/w*), or, 
equivalently, [y, n) -H- (— v, — n). The differential operator w* d w * flips the parity, so eq. (I4.16P 
should be odd under inversion. Since the two rational prefactors on the second line of eq. ( 14. 161) 
map into one another under inversion, we can infer that their coefficients must be relatecjj, 



fjo 



1 1 

w' w* 



-g (w,w 



gdw,w 



(4.17) 



where go(w,w*) = p(go(xo,xi)) and = p(g\(xo,xi)). It is easy to check that this 

identity is satisfied for low loop ordersj. 

Using these symmetry properties, we can write, 



w 



._d_ 
dw* 6 



MHVl 



2vri 



I LLA 



log(l — Mi) 1 + W* 



P 



W,W ) <H- 



1 1 

u>' w* 



(4.18) 



Turning to the right-hand side of eq. (I4.14p . we observe that the differential operator wd w acts 
on eq. (14. 3 P by removing the leading Xq and flipping the sign of the second term, 



3 ^,NMHV| 

ow 



2m 



LLA 



log(l - Ui) 1 



W* 



xz 



NMHV 



w w* 



(4.19) 



8 p does not generate any rational functions which might allow these terms to mix together. 

9 A general proof would be tantamount to showing that eq. (|4.2[) is symmetric under inversion. The latter 
seems to require another intricate cancellation of multiple zeta values. We postpone this investigation to future 
work. 
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Comparing eq. ( H~Tgj) and eq. fljTIgj) , we see that eq. ( dij) is satisfied if g- (^o,^i) = ^£ NMHV . 
To verify that this is true, we must extract go(x ,Xi) from -R^ hv |lla- To this end, collect all 
terms in the argument of p with at least one trailing xq and remove that xq. This procedure 
gives, 

i 00 / k ~ 2 n n2m-k+l \ 

— 2 \ ^ — 7TX — / 

_ p£ ^NMHV 

so we conclude that eq. (I4.14p is indeed satisfied. 



5 Collinear limit 

In the previous section, we proposed an all-orders formula for the MHV and NMHV remainder 
functions in MRK. The expressions are effectively functions of two variables, w and w*. The 
single-valuedness condition allows for these functions to be expressed in a compact way, but the 
result is still somewhat difficult to manipulate. 

In this section, we study a simpler kinematical configuration: the collinear corner of MRK 
phase space. To reach this configuration, we begin in multi-Regge kinematics and then take legs 
1 and 6 to be nearly collinear. In terms of the cross ratios Ui, this limit is 

1- Ul ,u 2 , « 3 ~0, Z = _^_ = 0(1) J y = ^^^Q } (5.1) 

1 — Ml 1 — Ui 

or, in terms of the (w,w*) variables, it is equivalent to, 

1 — u\ ~ , |w| ~ , w ~ w* . (5.2) 

As we approach the collinear limit, the remainder function can be expanded in powers of 
w, w*, and log | w | - The leading power-law behavior is proportional to (w + w*). Neglecting 
terms that are suppressed by further powers of |iu|, the result is effectively a function of a single 
variable, £ = r\ log \w\ = alog(l — u±) log \w\, and is simple enough to be computed explicitly, as 
we show in the following subsections. 



5.1 MHV 

In the MHV helicity configuration, the remainder function is symmetric under conjugation w <H> 
w*. It also vanishes in the strict collinear limit. These conditions suggest a convenient form for 
the expansion in the near-collinear limit, 



< HV | LL A,con. = . ™ J w + w*)J2v k+1 rr y (r}log\w\), (5.3) 

log(l - Ml) 
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for some functions r^ HV that are analytic in a neighborhood of the origin. We have neglected 
further power-suppressed terms, i.e. terms quadratic or higher in w or w*. The index k labels the 
degree to which r^ HV is subleading in log \ w\. For example, the leading logarithms are collected 
in r^ HV , the next-to-leading logarithms are collected in r]^ HV , etc. 

Starting from eq. (I4.2p . it is possible to obtain an explicit formula for r^ HV . To begin, we note 
that it is sufficient to restrict our attention to the terms proportional to w — the conjugation 
symmetry guarantees that they are equal to the terms proportional to w*. The main observation 
is that only a subset of terms in eq. (14. 2 p contribute to the power series expansion at order w. 
It turns out that the relevant subset is simply the set of SVHPLs with a single x\ in the weight 
vector. Roughly speaking, each additional x\ implies another integration by 1/(1 + w), which 
increases the leading power by one. 

The equivalent statement is not true for w*, i.e. SVHPLs with an arbitrary number of xi's 
contribute to the power series expansion at order w*. This asymmetry can be traced to the 
differences between the x and y alphabets: referring to eq. (I3.9p . the x alphabet indexes the 
HPLs with argument w and the y alphabet indexes the HPLs with argument w*. 

We are therefore led to consider the terms in eq. (14.21) with exactly one X\. Eq. (I4.4p shows 
that these terms may be obtained by dropping all X\S from X, 

pMHVi _ _ J |^?MHV 1 \ fr A \ 

R 6 lLLA,co,,.- log(1 _ ni) p(^ Z ~- 2 XlV). (5.4) 

Since no £ terms appear in SVHPLs with a single x±, it is straightforward to express them in 
terms of HPLS, 

n \ m \ 

jOt^t ,T m = / — tH'L H m n — i ~\- / — ~H n m-jH-. . (5.5) 
x xix / J .j x x% l X 1 X Q J / j 'J X%X!X X \ I 

j = J ' j=0 J ' 

Here we have simplified the notation by defining H m = H m (—w) and H m = H m (—w*). Next, we 
recall eq. (13.51) . in which we used the shuffle algebra to expose the explicit logarithms, 

H x n xixl n = J2 j^Sfji j J 11 "' J ■ ( 5 - 6 ) 

Finally, eqs. (13. 6 p and (13 .7p implies that the series expansions for small w have leading term, 

H x , Xl (-w) = -w + 0(w 2 ). (5.7) 

Combining eqs. (I5.4p -( 15T7|) and applying some hypergeometric function identities, we arrive 
at an explicit formula for r^ HV , 

k n 2k-n-m , _r,\2m+j-k-l 
«) = ^V + EE E L + ■ _ fc) , 3Mx™-^P^> k -^\0) I&yfr) . 
n=0 m=0 j=k—m 



k 

Z ' ■ ' ■ ' ■ ( III -f '/ — K I ! 

n=0 m=0 j=k—m v J ' 

(5.1 
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In this expression, the Ij are modified Bessel functions and the pi a ' b ^ are Jacobi polynomials, 
which can be defined for non-negative integers j by the generating function, 

oo , 

P t' b) ( z ) fj = T+b Vt 2 - 2tz + l) ~ a (l + t + y/t 2 - 2tz + l\ ~ (Vt 2 -2tz+l 



3=0 



(5.9) 

It is easy to extract the first few terms, 



r MHV (x) 



rf HV (x) 



i [1 - 7 (2>/5)] , 

-i/ 2 (2Vx), (5.10) 
(2V5)-l/ 4 (2Vi) . 



The leading term, r^ HV , corresponds to the double-leading-logarithmic approximation (DLLA) 
of ref. [3B], 

< HV |dlla = wr a (w + w*) [l-/ (2vVogM)] , (5.11) 

and is in agreement with the results of that reference. 

Only for k > 2 do £ values begin to appear in r^ HV . Moreover, modified Bessel functions 
with odd indices only appear in the ^-containing terms. To see this, notice that the £-free terms 
of eq. (15. 8 j) arise from the boundary of the sum with n = m = 0, in which case a = b = k — j in 
eq. (15. 9p . When a = b, Pj J ' b (0) = for odd j since eq. ( 15. 9ft reduces to a function of t 2 in this 
case. It follows that the £-free pieces of r^ HV have no modified Bessel functions with odd indices. 

Equations (15. 3p and ( 15. 8p provide an explicit formula for the six-point remainder function 
in the near-collinear limit of the LL approximation of MRK. If the sum in eq. (15. 3p converges 
sufficiently quickly, then it should be possible to evaluate the function numerically by truncating 
the sum at a finite value of k, /c max - A numerical analysis indicates that for \w\ < 1 and t] < 20, 
fc max ~ 100 is adequate to ensure convergence. 

The numerical analysis also indicates that -R^ hv |lla, coil, increases exponentially as a function 
of rj, and that the extent of this increase depends strongly on the value of log|w|. We find 
empirically that the rescaled function 

>MHVl _ ( V \ log(l- Ui) dMHV| 



^ 6 LLA, coll. - exp - == —— ■ — K 6 LLA, coll. {^-^ 

V ty- log \w\) 2m [W + W*) 

attains reasonable uniformity in the region < rj < 10 and —40 < log < 0. This particular 
rescaling carries no special significance, as alternatives are possible and may be more appropriate 
in different regions. In eq. (I5.12p we have also divided by the overall prefactor of eq. (I5.3P so that 
^g 11 ^ | lla, coil, is truly a function of the two variables rj and log \ w\. The results are displayed in 
Figure [TJ 
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o $> 



Figure 1: The MHV remainder function in the near-collinear limit of the LL approximation of 
MRK. It has been rescaled by an exponential damping factor. See eq. (I5.12p . 



5.2 NMHV 

A similar analysis can be performed for the NMHV helicity configuration. The situation is 
slightly more complicated in this case because the NMHV remainder function is not symmetric 
under conjugation w -R- w*. One consequence is that its expansion in the collinear limit requires 
two sequences of functions, which we choose to parameterize by r^ MHV and f^ MnY , 



dNMHVi 

-^6 |LLA, coll. 



2ni 



log(l - u Xt 



(w + w*)J2v k+2 ^ MUY {vlog\w\) 

oo 

+ ^£A~r HV (^og H ) 



fc=0 



(5.13) 



Contributions to the power series at order w arise from the first term of eq. (I4.3P (the second 
term has an overall factor of w*), and, as in the MHV case, only from the subset of SVHPLs 
with a single X\ in the weight vector. It is therefore possible to reuse eqs. f l5.5p -f l5T7|) and obtain 
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an explicit formula for the coefficient of w, r^ MHV . The result is, 



k n 2k—n—m , r\\2m+j-k 

= £E E ~ i ■ _ m, 3(n,m)x-'^-"/^y— (0) / J+1 (2VS) . 

n=0 m=0 j=k-m. \ J ) 

(5.14) 

1 



The first few terms are 



,NMHV/ 


NMHV 



=/i(2VS) 



L =/ 3 (2V5), (5.15) 



.r 



As previously mentioned, it is not so straightforward to extract the coefficient of w* in this 
way. We can instead make progress by exploiting the differential equation (I2.2ip . In terms of 
the functions r^ HV , r^ MHV , and f^ MHY , the equations read, 

d x rr\ X ) = 2 rr v (x) + d xr rrw r , 1 fi . 

^f fc NMHV (x) = 2r fc MHV (x) + 2rj™ v (x) • 

The first of these equations is automatically satisfied and confirms the consistency of eq. ( 15 .8p 
and eq. ( I5.14p . The second equation determines f^ MHV up to a constant of integration which can 
be determined by examining the n = —1 term of eq. (I2.19p . The solution is, 

k n 2k-n-m f ^m+j-k 
n=0 m=0 j=k-m 

(5.17) 

The first few terms are 



V^/i (2v^) 



-NMHV/ 
'0 



-NMHV, 



2 VXJ1 ^ VJ v ' (5.18) 
/ 1 (2 v / ^)-^/3(2v^) • 



2v^ 



Modified Bessel functions with even indices only appear in the ^-containing terms of r^ MHV and 
-NMHV rjkg eX pi ana tion of this fact is the same as in the MHV case, except that the parity is 
flipped due to the shifts of the indices of the modified Bessel functions in eq. (I5.14p and eq. (I5.17p . 
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5.3 The real part of the MHV remainder function in NLLA 

As described in Section the real part of the MHV remainder function in NLLA is related to 
its imaginary part in LLA. In the collinear limit, the relation (1 2 . 1 5 [) may be written as, 



Ke \ K 6 | NLLA, coll. J — , 7Z r I ~ V o — 77 log | W | J Kq | LLA, coll. 



2ni /l 2 d 1 1 
log(l - Mi) V2 77 dq f] 2 

log (1 - Mi) 

Since -R^ hv |lla vanishes like (w+w*) in the strict collinear limit, the quadratic term (-R^ hv |lla) 2 
in eq. (12.151) only contributes to further power-suppressed terms in the near-collinear limit and 
is therefore omitted from eq. f)5.19lj^l . We may write eq. (15.191) as, 

4vr 2 °° 

Re (# 6 mhv |nlla, con.) = Aw + w*) Vr/ fc+ V(r/log \w\) , (5.20) 

log (1 - Ui) ^ 

where, 

q k (x) = \x 5 , k + \ (k - x) rf^{x) + \ xd x rf^{x) . (5.21) 

The leading term, qo, corresponds to the real part of the next-to-double-leading-logarithmic 
approximation (NDLLA) of ref . jl6] . Our results agred^l with that reference and read, 

ReW HV |NDLL A ) = y- (™ + W *^ \-ri\og\w\I (2 ^ log \w\) + y/ rj log \w\h (2^gH 

loe: 1 — Mi L V / V 



(5.22) 



6 Conclusions 

In this article, we studied the six-point amplitude of planar M = 4 super- Yang-Mills theory in 
the leading- logarithmic approximation of multi-Regge kinematics. In this limit, the remainder 
function assumes a particularly simple form, which we exposed to all loop orders in terms of the 
single-valued harmonic polylogarithms introduced by Brown. The SVHPLs provide a natural 
basis of functions for the remainder function in MRK because the single-valuedness condition 
maps nicely onto a physical constraint imposed by unitarity. Previously, these functions had been 
used to calculate the remainder function in LLA through ten loops. In this work, we extended 
these results to all loop orders. 

In MRK, the tree amplitudes in the MHV and NMHV helicity configurations are identical. 
This observation motivates the definition of an NMHV remainder function in analogy with the 

10 As a consequence, eq. f)5. 19|) does not depend on the conventions used to define R, i.e. the equation is equally 
valid if R is replaced by exp(i?). 

n The agreement requires a few typos to be corrected in eq. (A. 16) of ref. [46] , 
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MHV case. We examined both remainder functions in this article, and proposed all-order formulas 
for each case. In fact, these formulas are related: as described in ref. [39} , the two remainder 
functions are linked by a simple differential equation. We employed this differential equation to 
verify the consistency of our results. 

We also investigated the behavior of our formulas in the near-collinear limit of MRK. The 
additional large logarithms that arise in this limit impose a hierarchical organization of the 
resulting expansions. We derived explicit all-orders expressions for the terms of this logarithmic 
expansion. The results are given in terms of modified Bessel functions. 

We did not provide a proof of the all-orders result, but we verified that it agrees through 
14 loops with an integral formula of Lipatov and Prygarin. The agreement of these formulas 
at 12 loops and beyond requires an intricate cancellation of multiple zeta values. It would be 
interesting to understand the mechanism of this cancellation. There are several other potential 
directions for future research. For example, in refs. jS3H53], Alday, Gaiotto, Maldacena, Sever, 
and Vieira performed an OPE analysis of hexagonal Wilson loops which in principle should 
provide additional cross-checks of our results. It should also be possible to study the all-orders 
formula as a function of the coupling and, in particular, to examine its strong-coupling expansion. 
We have begun this study in the collinear limit and presented our initial results in Figure [TJ A 
first attempt to compare the six-point remainder function in MRK at strong and weak coupling 
was made by Bartels, Kotanski, and Schomerus [32]. Further analysis of our all-orders formula 
should allow for an important comparison with this string-theoretic calculation. 
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